LUDMILA MARIA RATTIS TEIXEIRA

Exerccio 5:

Usando o conjunto de dados Conjunto de Dados: Dados de Biomassa de rvores de Eucalyptus Saligna, construa os seguintes grficos: 

##Criando e editando um grfico

Um grfico de disperso entre dap e ht editando:

1-)Legenda em ambos os eixos

2-)Marcaes (ticks) para dentro 

3-)Com apenas dois eixos (formato L)

4-)Ttulo do grfico

R:
  plot(esaligna$dap~esaligna$ht, xlab="altura", ylab="DAP", bty="l", tcl=0.3,
  title(main="Altura de E.saligna em funo do DAP"))


5-)Tamanho de fonte

R:
  par(cex.axis=1.5)
  par(cex.lab=2)


###################################################################

##Construindo dois grficos

Dois grficos lado a lado entre dap e talhao:

1-)Construir um grfico de desenho de caixa (boxplot) entre as variveis
dap e classe e ao lado um plot com apenas a mdia e desvio-
padro. E insira tambm uma letra para dizer qual  o grfico a e
qual  o b (tanto faz, quem  um e quem  outro). 

##Problema no enunciado

esaligna$classe = as.factor(esaligna$classe)
classe = sort(unique(esaligna$classe))
classe = as.data.frame(classe)
media = tapply(esaligna$dap,esaligna$classe,mean)
sd = tapply(esaligna$dap,esaligna$classe,sd)
classe$media = as.vector(media)
classe$sd = as.vector(sd)
classe
par(mfrow=c(1,2))
boxplot(esaligna$dap~esaligna$classe, main="Boxplot DAP X Classes")
text(15,15,"a",cex=1.8)
plot(classe$media~classe$sd, ylab="Media das Classes", xlab="Desvio das Classes", main="Plot Desvio X Classes", bty="L")
text(15,15,"b",cex=1.8)
par(mfrow=c(1,1))

Dica: vocs tero que calcular a mdia e os devio padro do dap das
rvores em cada classe. Depois crie uma matriz com estes valores e
crie o plot destes valores. 
###################################################################

##Adivinhando o cdigo

1-)Leia o arquivo exercicio3 e usando os dados x1 e y1 e x2 e y2
tente reproduzir esta figura: 

R:

##Grfico 1

desc=read.table("exercicio3.csv",as.is=T,header=T,sep=",")
par(mfrow=c(1,2))
str(desc)
plot(desc$y1~desc$x1,xlab="Log(PatchSize)ha",ylab="Euclidean distances",xaxp=c(0.5,2.5,5),yaxp=c(0,3,8),bty="l",pch=17)
model=lm(desc$y1~desc$x1)
abline(model,lwd=1.5)
text(3,3,"a",cex=3.8)

##Grfico 2

legenda.x = as.vector(c("Small","Medium \n Edge","Medium \n Interior","Large \n Edge","Large \n Interior","Control"))
par(bty="l")
par(cex.lab=1)
par(cex.axis=1)
par(mar=c(5,4,5,1)) 
boxplot(desc$x2~desc$y2,names=legenda.x, cex.axis=0.7, ylim=c(0,3))
text(7,3,"b",cex=1,5)
text(1,3,"*",cex=2)
text(2,3,"*",cex=2)
text(3,3,"**",cex=2)
text(4,3,"*",cex=2)
text(5,3,"***",cex=2)
par(mfrow=c(1,1)



 


